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We perform a systematic study of the possible loosely bound states composed of two charmed 
baryons or a charmed baryon and an anti-charmed baryon within the framework of the one boson 
exchange (OBE) model. We consider not only the tt exchange but also the t], p, lo, cf) and o exchanges. 
The S — D mixing effects for the spin-triplets are also taken into account. With the derived effective 
potentials, we calculate the binding energies and root-mean-square (RMS) radii for the systems 
AcAc(Ac), HcHc(Hc), EcEc(Ec), H^H^H^) and Q,cO,c{^c)- Our numerical results indicate that: (1) 
the H-dibaryon-Iike state AcAc does not exist; (2) there may exist four loosely bound deuteron-like 
states HcHc and HcHj, with small binding energies and large RMS radii. . 



PACS numbers: 12.39.Pn, 14.20.-c, 12.40.Yx 
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Many so-called "XYZ" charmonium-like states such as X(3872), X(4350) and F(3940) have been observed by Belle, 
f~i ' CDF, DO and BaBar collaborations [1-4] during the past fe-w years. Despite the similar production mechanism, some 
Oh- of these structures do not easily fit into the conventional charmonium spectrum, -which implies other interpretations 

such as hybrid mesons, heavy meson molecular states etc. might be responsible for these ne-w states [5] [6]. 
^ , A natural idea is that some of the "XYZ" states near two heavy meson threshold may be bound states of a pair 
' of heavy meson and anti-heavy meson. Actually, Rujula et al. applied this idea to explain ■!/'(4040) as a P--wave 
D*D* bound resonance in the 1970s [7]. Tornqvist performed an intensive study of the possible deuteron-like t-wo- 
' charm-meson bound states -with the one-pion-exchange (OPE) potential model in Ref. [8]. Recently, motivated by the 
J> controversy over the nature of X(3872) and Z(4430), some authors proposed X{3872) might be a DD* bound state [9- 
f — 13]. Our group have studied the possible molecular structures composed of a pair of heavy mesons in the frame-work 
of the One-Boson- Exchange (OBE) model systematically [14, 15]. There are also many interesting investigations of 
other hadron clusters [16-21]. 

' The boson exchange models are very successful to describe nuclear force [22-24] . Especially the deuteron is a loosely 
, bound state of proton and neutron, -which may be regarded as a hadronic molecular state. One may -wonder -whether 
' a pair of heavy baryons can form a deuteron-like bound state through the light meson exchange mechanism. On the 
j I other hand, the large masses of the heavy baryons reduce the kinetic of the systems, -which makes it easier to form 
" . l bound states. Such a system is approximately non-relativistic. Therefore, it is very interesting to study -whether the 
I> OBE interactions are strong enough to bind the t-wo heavy baryons (dibaryon) or a heavy baryon and an anti-baryon 
(baryonium) . 

A heavy charmed baryon contains a charm quark and t-wo light quarks. The f-wo light quarks form a diquark. 
, Heavy charmed baryons can be categorized by the flavor -wave function of the diquark, -which form a symmetric 6 or 
an antisymmetric 3 representation. For the ground heavy baryon, the spin of the diquark is either or 1, and the spin 
of the baryon is either 1/2 or 3/2. The product of the diquark flavor and spin -wave functions of the ground charmed 
baryon must be symmetric and correlate with each other. Thus the spin of the sextet diquark is 1 while the spin of 
the anti-triplet diquark is 0. 

The ground charmed baryons are grouped into one antitrpilet with spin- 1/2 and two sextets with spin- 1/2 and 
spin-3/2 respectively. These multiplets are usually denoted as B^, Bq and Bq in literature [25]. In the present work, 
we study the charmed dibaryon and baryonium systems, i.e. AcAc(Ac), 5cSc(Sc), ScSc(Sc), '^'c'^'d'^'c) ^^'^ ^c^ci^c)- 
Other configurations will be explored in a future work. We first derive the effective potentials of these systems. Then 
we calculate the binding energies and root-mean-square (RMS) radii to determine which system might be a loosely 
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bound molecular state. 

This work is organized as follows. We present the formalism in section II. In section III, we discuss the extraction 
of the coupling constants between the heavy baryons and light mesons and give the numerical results in Section IV. 
The last section is a brief summary. Some useful formula and figures are listed in appendix. 

II. FORMALISM 

In this section we will construct the wave functions and derive the effective potentials. 

A. Wave Functions 

As illustrated in Fig. 1, the states A+, S+ and S° belong to the antitriplet Bg while S^+, E!^° and fi^ 

are in sextet Bg- Among them, A+ and 17° arc isoscalars; {S+, S°} and E'^} are isospin spinnors; {S^"*", S+, E°} 
is an isovector. We denote these states Ag, S^, Sg, and ^c- 




(a) antitriplet (b) sextet 

FIG. 1: The antitriplet and sextet. Here the brackets and parentheses represent antisymmetrization and symmetrization of the 
light quarks respectively. 

The wave function of a dibaryon is the product of its isospin, spatial and spin wave functions, 

*L'/;'''^"-*L«*L«'^L- (1) 

We consider the isospin function VP^^ first. The isospin of Ac is 0, so AcAc has isospin 7 = and ^'i^Ae = A+A+, 
which is symmetric. For Sc5c, the isospin is 7 = or 1, and their corresponding wave functions are antisymmetric 
and symmetric respectively. Ec^c has isospin 0, 1 or 2. Their flavor wave functions can be constructed using Clebsch- 
Gordan coefficients. S^S^ is the same as ScSc. The isospin of the ^c^c is 0. Because strong interactions conserve 
isospin symmetry, the effective potentials do not depend on the third components of the isospin. For example, it is 
adequate to take the isospin function with 73 = 1 when we derive the effective potential for 'i'lr- , though the 

wave function ;^(S+S° + S°5+) indeed gives the same result. In the following, we show the relevant isospin functions 
used in our calculation, 
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Wc arc mainly interested in the ground states of dibaryons and baryonia where the spatial wave functions of these 
states are symmetric. The tensor force in the effective potentials mixes the S and D waves. Thus a physical ground 
state is actually a superposition of the S and D waves. This mixture fortunately does not affect the symmetries of 
the spatial wave functions. As a mater of fact, for a dibaryon with a specific total spin J, we must add the spins of its 
components to form S first and then couple S and the relative orbit angular momentum L together to get J = L + S. 
This L — S coupling scheme leads to six S and D wave states: ^5*0, '^5'i, ^D2, "^fi, ^D^ and ^D^- But the tenser force 
only mixes states with the same S and J. In our case we must deal with the ^S\-^Di mixing. After stripping off the 
isospin function, the mixed wave function is 

\i^)=Rs{r)\^S{)+RD{T)\^D^), (7) 

which will lead to coupled channel Schrodinger equations for the radial functions Rs{r) and Roir). In short, for the 
spatial wave functions, we will discuss the ground states in ^Sq and ^S\, and the latter mixes with ^Di. 

Finally, wc point out that the / and S of states in Eq. (1) can not be combined arbitrarily because the generalized 
identity principle constricts the wave functions to be antisymmetric. It turns out that the survived compositions are 
*aI1' ^sJe' *HchI' *eI1'*eIec'*sI1'*h^h\' *h^hm and ^^^l]^. For baryonia, there is no constraint on the wave 
functions. So we need take into account more states. The wave functions of baryonia can be constructed in a similar 
way. However, we can use the so-called "G-Parity rule" to derive the effective potentials for baryonia directly from 
the corresponding potentials for dibaryons, and it is no need discussing them here now. 



B. Lagrangians 

We introduce notations 

Ac = A+ Se=|-^(-Sj+ + E°),-J=(-E++-SO),E+|, , $7^ = 0° (8) 

to represent the corresponding baryon fields. The long range interactions are provided by the tt and r] meson exchanges: 

-Ctt = 5'7rB<,Bc2c«75TSc • TT + ff^ScS, {~i)'^cil5 X Sc ' + 9TTE'^E'^%i'y5T% ■ TT (9) 
+9r,^c^c'^c ■ hs^cV + griE'^E'^KhsKV + griQcQc^chb^cV, (10) 

where Qt^ScBc^ STrEcSc) etc. are the coupling constants, r = {Ti,r2,T3} are the Pauli matrices, and tt = 

{;^(7r+ + 7r~), -^{t^^ — tt"), tt'^} are the tt fields. The vector meson exchange Lagrangians read 

+5.HiHiS;7^TS; . p^ + :teks;a,.TS; • d^p'^ (11) 
+g^s'^E'M^,K^' + ^S-c^^^-od^^" (12) 
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+94,n^a.^clt.^cr + ^-^^n.a^M^r, (13) 
with p — {-^(p^ + P )i ^'"^ ^ P ): The a exchange Lagrangian is 

+5^s^Hi5c2c'^ + 5^n,n.f^cf^cCT. (14) 
There are thirty-three unknown couphng constants in the above Lagrangains, which will be determined in Sec. III. 



C. Effective Potentials 



To obtain the effective potentials, we calculate the T matrices of the scattering processes such as Fig. 2 in momentum 
space. Expanding the T matrices with external momenta to the leading order, one gets [26] 

Vir) = J^I d\e-^Q-^T{Q)F{Q)\ (15) 

where J-{Q) is the form factor, with which the divergency in the above integral is controlled, and the non-point-like 
hadronic structures attached to each vertex are roughly taken into account. Here we choose the monopole form factor 

with Q — {Qo, Q} and the cutoff A. 



At 



At 



A+ A+ A+ K 

FIG. 2: Scattering processes of AcAc — > AcAc and AcAc AcAc, Qs are the transformed four momenta. 



Generally speaking, a potential derived from the scattering T matrix consists of the central term, spin-spin inter- 
action term, orbit-spin interaction term and tenser force term, i.e., 

V{r) = Vcir) + Vss{r)(Ti ■ ^2 + VLs{r)L ■ S + VT{r)Si2{r), (17) 

where Si2{r) is the tensor force operator, Si2{r) = 3((Ti • f){(T2 ■ r) — eri ■ (T2. The effective potential of a specific 
channel, for example AcAc AcAc shown in Fig. 2, may contain contributions from the pseudoscalar, vector and 
scalar meson exchanges. We need work them out one by one and add them. The potentials with the stripped isospin 
factors from the pseudoscalar, vector and scalar (cr here) meson exchange are 

y"(r;a,/i) = ^^"^(r; a, /i)^i • 0-2 + l/^(r; a, /i)5i2(r), 

V'ir; /?, h) = vSir; /?, h) + V^si^l /?, h)a, ■ ^2 + Vlsir; /3, h)L ■ S + ^^(r; /3, /i)5i2(r), 

V'{r; a, h) = F^(r; a, h) + V£sir; a, h)L ■ S, (18) 



5 



where a = Tr,r], = p,u),(f) and 



Vs^sir;a,h) = - 



9lhh ml 
47r 12m| 



Hi{K,ma,r), 



1 



3hh 



Ho{A,mi3,r) - {g^^^f^ + 4g,shhfi3hh)-^Hi{A,mp,r) 



rrir. 



-^{9fihh + fphh) -^HiiA, mi3,r), 



Vls{r\P,h) = - — {3gjf^f^ + Ag^hhf0hh)^H2{A,m0,r 



VS{r;(3,h) 
VS{r;a,h) 
V£sir;a,h) 



1 

An 



-jzi9/3hh + ffshh) j^H3{A,mf},r), 



-m. 



9ahh 

An 



HoiA,m^,r) + —^Hi{A,m^,r) 



(19) 



The definitions of functions Hq, Hi, H2 and are given in the appendix. Prom Eq. (18), one can see the tensor 

force terms and spin-spin terms arc from tlic pscudoscalar and vector meson exchanges while the central and obit-spin 
terms are from the vector and scalar meson exchanges. Finally the effective potential of the state hh is 

Vhh{r) = Y,eiV%r;a,h) + Y,ClV\r;^,h) +e'',V'^{r;a,h) 

a 

Y,ClV^{r;P,h)+ClV§{r;a,h) \ + J ^^^^^^(r; a, /i) + ^C^F|s(r; /3, /i) \<J^ a2 



+ { Y.C'ffVUr; 13, h) + C^VEsir; P,h)\L-S- 



y^VSir; a, h) + ^C^l4(r; /3, h) ) S^^ir), 



(20) 



where C^, and are the isospin factors, which are listed in Table I. 
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TABLE I: Isospin factors. The values in brackets for baryonia are derived by the "G-Parity rule". 



Given the effective potential 14^, the potential for hh, V^^, can be obtained using the "G-Parity rule", which states 
that the amplitude (or the effective potential) of the process AA — >• AA with one light meson exchange is related 
to that of the process AA AA by multiplying tlie latter by a factor (— )^'^, where (— )^° is the G-Parity of the 
exchanged light meson [27]. The expression of V^^ is the same as Eq. (20) but with V°'{r;a,h), V^{r;p,h) and 
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y^(r; cr, h) replaced by V°'{r; a, h), V^{r; /3, h) and V'^{r; a, h) respectively. 

V''{r;a,h) = (-)^°MF''(r; a, /i), 
V\r;p-h) = {-y'^Wv'>{r;p,h), 

V%r;a,h) = (-)^<^My'=(r; tr, /i). (21) 

For example, 

y«(r;a;,Ae) = (-l)F«(r;u;, A^), (22) 

since the G-Parity of oj is negative. In other words, we can still use the right hand side of Eq. (20) to calculate 
but with the redefined isospin factors 

C«^(-)^oMC«, ^ (-)^°[«C^, C^^(-)^°WC^, (23) 

which are listed in Table I too. 

The treatments of operators ai ■ (T2, L ■ S and 5i2(r) are straightforward. For ^Sq, 

(71 •0-2 = -3, L-S = 0, 5i2(r)=0, (24) 

which lead to single channel Shrodinger equations. But for ^^i, because of mixing with ^Di, the above operators 
should be represented in the {pSi), |^-Di)} space, i.e., 

.... ^ (;;),..«= (2 _°3),s„,.).(^ 5)^ ,25, 

These representations lead to the coupled channel Shrodinger equations. 

III. COUPLING CONSTANTS 

It is difiicult to extract the coupling constants in the Lagrangians experimentally. We may estimate them using the 
well-known nucleon-meson coupling constants as inputs with the help of the quark model. The details of this method 
are provided in Ref. [28] . The one-boson exchange Lagrangian at the quark level is 

+gmq [uilburj + di-fsdrj - 2sij5sr]) 

+9pqq {uit^Up"'' - d-ff.dp"'') 

+9<^qq {u'Jij.UUJ^' + d'yi^duj^) + g4,qqS^^S(j)'^ 

+gaqq [uucF + McF + sso) H , (26) 

where g-^qq, gnqq, ■ ■ g^qq are the coupling constants of the light mesons and quarks. The vector meson terms in 

this Lagrangian do not contain the anomalous magnetic moment part because the constituent quarks are treated as 
point-like particles. At the hadronic level, for instance, the nucleon-nuclcon-meson interaction Lagrangian reads 

jCnn = gwNNNi-f^rN ■ tt + gnNNNij^Nr] 

+gpNNNl,.TN ■ + IfUNa^^rN ■ d^p^ 



fujNN 

2m N 

g^NNNNa, (27) 



where g^NN, g-qNN, ■ ■ ■, gaNN are the coupling constants. We calculate the matrix elements for a specific process both 
at quark and hadronic levels and then match them. In this way, we get relations between the two sets of coupling 
constants, 

5 niN tun 
3^ 
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_ q 9u>NN + fojNN _ 9uiqq 

9u}NN — oduiggj — ) 

mjv niq 

_ 9pNN + fpNN _ 5 Qpqq 

9pNN - 9pqq, — - o > 

THj^ O TTlq 

9aNN = Sgcrqq. (28) 

From these relations, we can see that q^nn and /^nn are not independent. So are QpNN and fpNN- The constituent 
quark mass is about one third of the nucleon mass. Thus we have fujNN ~ and fpNN ~ ^9pNN- 

With the same prescription, we can obtain similar relations for heavy charmed baryons which are collected in the 
appendix. Substituting the coupling constants at the quark level with those from Eq. (28), we have 

4 ms^ 2 mn^ 

5^=0=0=0) SttEcEc = cflwiVAr , 9ttE'S' = -^g%NN -, (29) 

5 rriN " " 5 mjv 
n n 4 ms, 

2 '^H^ 8 mn^ 

9riS'E' — —i:9r)NN , g-nQ^Q^ = --^griNN , (,oUj 

o rriN 6 TTlN 

2 2 2 

gaAcAc = -^gaNN, S'ffHeHe = -^gcrNN, So-EeEe = gScrATiV, 

2 2 

SctHJ^H^ = -^gaNN, gcrncQc = -^g^NN, (31) 

2 , 2 

ffwAcAc — -^gujNN, JwAcAo — —-^g^NN, 

1 , 1 

gwScSa — -^guNN, JuiSaB., = -gfftjATJV, 

2 . 2 / rUE, ^ 
5wEcEc = i;guNN, JwEcEc — i;guNN ^ J- 

6 6 \ mjv 

1 1 / 

5''«'=o=i = gflc^ATJV, /u-H^H^ = gSftuJVAT (2—=^ - 1 ) , (32) 

gp-n^T,^ = "^gpNN, fpS^T,^ = 3 (SpJVAT + /piVJv) (2^^ - 1 ) , 

V mjv / 

1 / Tfl"' \ 

gpE'^E'^ = gpNN, fpE'^E'^ = -^{gpNN + fpNN) - Ij , (33) 

g<l>EcEc = V^gpNN, UEcEc = {gpNN + /pJVAr) , 



g<l>E'^E'^ = V^gpNN, f<pE'^E'^ = -3- [gpNN + fpNNJ I 2 ^ — 1 



g<t>Qa^c 



= 2V2gpNN, UqMc = -T~i9pNN + fpNN) i 2—^ - 1 ) ) (34) 



where wc have used tun ~ Sm^. The couplings of <j) and heavy charmed baryons can not be derived directly from the 
results for nucleons. So in the right hand side of Eq. (34), we use the couplings of p and nucleons. 

The above formula relate the unknown coupling constants for heavy charmed baryons to g^NN, griNN, etc. which 
can be determined by fitting to experimental data. We choose the values g^NN = 13.07, grjNN = 2.242, gaNN = 8.46, 
guiNN = 15.85, fujNN / g^NN = 0, gpNN = 3.25 and fpNN/gpNN = 6.1 from Refs. [22, 23, 29] as inputs. In Table II, we 
list the numerical results of the coupling constants of the heavy charmed baryons and light mesons. One notices that 
the vector meson couplings for Sc^c and AgAc have opposite signs. They almost cancel out and do not contribute 
to the tensor terms for spin-triplets. Thus in the following numerical analysis, we omit the tensor forces of the 
spin-triplets in the ScSg and AgAc systems. 
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5.64 
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5.64 


5.64 




5.64 


UJ 


10.57 -10.57 


5.28 


-5.28 


10.57 44.67 


5.28 


23.72 




P 




3.25 


-4.62 


6.50 39.01 


3.25 


20.72 




</> 




4.60 


-6.53 




4.60 


29.30 


9.19 61.94 



TABLE II: Numerical results of the coupling constants. The coupling constants with the 4> exchange are deduced from Qpnn- 

IV. NUMERICAL RESULTS 

With the effective potentials and the coupling constants derived in the previous sections, one can calculate the 
binding energies and root-mean-square (RMS) radii for every possible molecular state numerically. Here we adopt 
the program FESSDE which is a FORTRAN routine to solve problems of multi-channel coupled ordinary differential 
equations [30]. Besides the coupling constants in Table II, we also need heavy charmed baryon masses listed in 
Table III as inputs. The typical value of this cutoff parameter for the deuteron is 1.2 ~ 1.5 GeV [22]. In our case, the 
cutoff parameter A is taken in the region 0.80 ^ 2.00 GeV. Such a region is broad and reasonable enough to give us 
a clear picture of the possibility of the heavy baryon molecules. 



baryon 


mass(MeV) baryon 


mass(MeV) 


meson 


mass(MeV) 


meson 


mass(MeV) 


At 


2286.5 


Ec 


2455 


7r± 


139.6 


P 


775.5 


■— 'C 


2467.8 


' — 'c 


2575.6 




135.0 


UJ 


782.7 


■— 'c 


2470.9 


■~I0 


2577.9 


V 


547.9 


4> 


1019.5 




2695.2 






a 


600 







TABLE III: Masses of heavy baryons and light mesons [31]. We use ms^ = 2469.3 MeV, to-h^ = 2576.7 MeV and rrin = 
138.1 MeV as numerical analysis inputs. 



A. AcAc and HcHc systems 

The total effective potential of AcAc arises from the a and w exchanges. We plot it with A = 0.9 GeV in Fig. 3 (a), 
from which we can see that the w exchange is repulsive while the a exchange is attractive. Because of the cancellation, 
the total potential is too shallow to bind two AgS. In fact, we fail to find any bound solutions of ^'Iv']^ even if one 
takes the deepest potential with A = 0.9 GeV. In other words, the loosely bound AgAg molecular state does not exist, 
which is the heavy analogue of the famous H dibaryon [32-35] to some extent. 

For the AcAc system as shown in Fig. 3 (b), both a and u exchanges are attractive. They enhance each other and 
lead to a very strong total interaction. From our results listed in Table IV, the binding energies of the AgAc system 
could be rather large. For example, when we increase the cutoff to A = 1.10 GeV, the corresponding binding energy 
is 142.19 MeV. The binding energies and RMS radii of this system are very sensitive to the cutoff, which seems to be 
a general feature of the systems composed of one hadron and anti-hadron. 

2° and S+ contain the s quark and their isospin is I = 1/2. Besides the a and w meson exchanges, the </> and p 
exchanges also contribute to the potentials for the ScSc(Sc) systems. Figs. 3 (c) and (d) illustrate the total potentials 
and the contributions from the light meson exchanges for ^B'^' and . For ^[^'^ , the attraction arises from the 

cr exchange. Because of the repulsion provided by the 0, p and oj exchange in short range, the total potential has a 
shallow well at r « 0.2 fm. However, the exchange almost does not contribute to the potential of ^^'S and the p 
exchange is attractive which cancels the repulsion of the a exchange. The total potential is about two times deeper 
than the total potential of . 

In Table IV, one notices that the binding energy is only hundreds of keV for ^'i^'ij when the cutoff varies from 
1.01 GeV to 1.20 GeV. Moreover the RMS radius of this bound state is very large. So the state ^'3'= is very 
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r|GcV 'l rlGeV'l r |GeV '| 



(a) yi^'^i with A = 0.9 GeV. (b) vj^y^ with A = 0.9 GeV. (c) VP'sj^ with A = 1.1 GeV. 




02468 10 02468 10 02468 10 

r[GeV"'] r[GeV-'] r[GeV-'] 



(d) l4°;|> with A = 1.0 GeV. (e) fJ"'^' with A = 1.0 GeV. (f) vl^J?[ with A = 1.0 GeV. 

FIG. 3: The potentials of '^AcAc, ^AoAoj ^ScSc and ^fg^g^. The spin-triplets have no S — D mixing because of the cancellations 
of the coupling constants. 
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2.53 
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1.00 


7.41 
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4.48 
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0.96 




1.10 


0.35 


4.62 
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0.78 




1.20 


0.18 


5.40 




0.87 


1.48 


2.72 




0.90 


1.24 


2.92 


^[0.1(3)] 


0.90 


4.12 


1.78 


^[1,1(3)1 


1.00 


10.33 


1.25 




1.00 


28.94 


0.86 




1.10 


31.80 


0.83 




1.10 


82.86 


0.60 




1.20 


66.19 


0.64 



TABLE IV: Numerical results of the systems AcAc, AcAc, — cHc and ^cHc, where "— " means no bound state solutions. After 
neglecting the tensor force terms, the results of spin-triplets are the same as those of spin-singlets. 



loosely bound if it really exists. The W~'J bound state may also exist. Its binding energy and RMS radius are 
2.53 - 36.59 MeV and 2.17 - 0.78 fm resp'ee-tively with A = 0.95 - 1.20 GeV. 

As for the Sc5c systems, the potentials are very deep. The contribution from the (p exchange is negligible too, as 
shown in Fig. 3 (e) and (f). We find four bound state solutions for these systems: , ^l/^'i , ^^'^ and ^B'g • 

Among them, the numerical results of and are almost the same as those of "^-^^ and respectively. 

The binding energies and the RMS radii of these states are shown in Table IV. We can see that the binding energy of 
\E'!?'i' varies from 1.48 MeV to 82.86 MeV whereas the RMS radius reduces from 2.72 fm to 0.60 fm when the cutoff 

is below 1.10 GeV. The situation of is similar to that of ^-''^ qualitatively. They may exist. But the binding 

energies appear a little large and the RMS radii too small when one takes A above 1.10 GeV. 



B. ScEc, H^H^ and ficfic systems 

For the SgSg system, all the tt, 77, c, lu and p exchanges contribute to the total potential. We give the variation of 
the potentials with r in Figs. 4 (a) and (b). For ^'^'^ > the potential of the uj exchange and p exchange almost cancel 
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out, and the 77 exchange gives very smaU contribution. So the total potential of this state mainly comes from the 
TT and a exchanges which account for the long and medium range attraction respectively. There may exist a bound 
state ^^^^ji^, see Table V. 

But for the other spin-singlet, > exchange provides only as small as 0.2 GeV attraction while the w and 

[2 11 

p exchanges give strong repulsions in short range r < 0.6 fm. We have not found any bound solutions for ^'^'ji as 

[1 31 

shown in Tabic V. For the spin-triplet state ^'^'s , there exist boimd state solutions with binding energies between 
0.11 MeV and 31.35 MeV when the cutoff lies between 1.05 GeV and 1.80 GeV. This state is the mixture of ^Si and 
^Di due to the tensor force in the potential. Prom Table V, one can see the S wave percentage is more than 90%. 




(a) V^"'^' with A = 1.10 GeV. 



(d) V'^'^' with A = 1.0 GeV. 
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(f) V^,'l\ with A = 1.70 GeV. 
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(i) with A = 1.0 GeV. 



J-'"/ 



(j) Vj^°'il^ with A = 1.0 GeV. 



FIG. 4: Potentials of ^e^e,, *e,e,, *s^a',, ^'a^s^, ^ficHc and *n^n^ 
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^[0.1] 



[1.1] 



[2,1] 
EcEc 



A (GcV) E (MeV) r™s (fm) 



1.07 


1.80 


2.32 


1.08 


3.10 


1.88 


1.10 


6.55 


1.44 


1.20 


42.95 


0.78 


1.25 


75.75 


0.65 



[0,3] 



1,3] 



:2,3] 



A (GeV) E (MeV) rrms (fm) Ps : Pd (%) 



1.05 
1.47 
1.50 
1.80 



0.11 
2.03 
2.52 
31.35 



5.94 
2.48 
2.27 
0.76 



98.11 1.89 

94.21 5.79 

93.79 6.21 

91.41 8.59 



TABLE V: 
solutions. 



Numerical results of the EcEc system, where the symbol "x" means this state is forbidden and 



means no 





One Boson Exchange 


One Pion Exchange 




A(ljcV) 


(MeV J 


rrms(im) 


Ps ■ 


Pd (% ) 


A(GeV) 


E(MeV) 


rrms(im) 


Ps 


Pd (%) 




0.97 


0.86 


3.76 
















^[(),1J 

EcEc 


0.98 


3.03 


2.21 


















1.00 
1.05 


18.43 
175.56 


1.01 
0.41 


















0.93 


1.04 


3.50 


81.20 


18.80 


0.80 


17.54 


1.20 


82.93 


17.07 


^[0,3] 
EoEc 


0.94 


2.55 


2.57 


75.27 


24.73 


0.85 


26.33 


1.04 


81.66 


18.34 




1.00 


28.16 


1.29 


58.07 


41.93 


0.90 


37.48 


0.92 


80.57 


19.42 




1.05 


78.48 


0.99 


50.56 


49.44 


1.05 


87.94 


0.68 


78.03 


21.97 




0.93 


0.75 


3.77 
















*e'j: 


0.94 


2.54 


2.27 


















0.98 
1.00 


32.28 
66.97 


0.80 
0.60 


















0.80 


3.71 


1.91 


94.73 


5.27 


0.97 


1.04 


3.14 


93.68 


6.32 


^[1,3] 
EcEc 


0.81 


5.18 


1.69 


94.38 


5.62 


1.02 


2.51 


2.18 


91.58 


8.42 




0.90 


40.35 


0.86 


90.12 


9.88 


1.10 


6.44 


1.51 


89.04 


10.96 




1.00 


143.46 


0.62 


76.86 


23.14 


1.30 


27.27 


0.88 


84.89 


15.11 


*E 


0.80 


24.87 


0.85 






0.75 


2.49 


1.98 






0.85 


49.30 


0.67 






0.80 


5.95 


1.38 








0.90 


90.04 


0.55 






0.90 


18.30 


0.88 








0.95 


149.66 


0.46 






1.10 


72.23 


0.51 






^[2,3] 
EoEc 


0.90 


1.44 


2.93 


96.92 


3.08 












1.00 


14.99 


1.21 


95.43 


4.57 














1.10 
1.20 


41.81 
77.28 


0.86 
0.71 


95.11 
94.72 


4.89 
5.28 













TABLE VI: Numerical results of the Ec-Ec system. Results from the OBE and OPE alone are compared. 



There exist bound state solutions for all six states of the Sc^c system. The potentials of the three spin-singlets are 
plotted in Figs. 4 (c)-(e). The attraction that binds the baryonium mainly comes from the p and w exchanges. These 
contributions are of relatively short range at region r < 0.6 fm. One may wonder whether the annihilation of the 
heavy baryon and anti-baryon might play a role here. Thus the numerical results for ScSc with strong short-range 
attractions should be taken with caution. This feature differs from the dibaryon systems greatly. 

In Table VI, for comparison, we also present the numerical results with the tt exchange only. It's very interesting 
to investigate whether the long-range one-pion-exchange potential (OPE) alone is strong enough to bind the baryonia 

and form loosely bound molecular states. There do not exist bound states solutions for V&^'ij and since the tt 

exchange is repulsive. In contrast, the attractions from the tt exchange are strong enough to form baryonium bound 
states for , f^i'^ and ^^'^ . We notice that the 5 - mixing effect for the spin-triplets mentioned above is 

stronger than that for the S^Sc system. 
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A (GeV) E (MeV) Vrms (fm) 



A (GeV) E (MeV) 



(fm) Ps : Pd (%) 



[0,3] 



0.95 


1.22 


3.03 


97.88 


2.12 


0.98 


2.44 


2.29 


97.45 


2.55 


1.00 


3.41 


2.01 


97.26 


2.74 


1.20 


15.43 


1.16 


96.74 


3.26 


1.30 


21.50 


1.03 


96.83 


3.17 



[1.1] 



1.50 
1.65 

1.70 

1.80 
1.90 



0.18 
1.24 

1.83 

3.42 
5.58 



5.52 
3.08 

2.64 

2.08 
1.74 



:i,3] 



TABLE VII: Numerical results of the 2^2^ system. 







One Boson Exchanj 


^es 




One Pion Exchanges 




A (GeV) 


E (MeV) 


rrma (fm) 


Ps : 


Pd{%) 


A (GeV) 


E (MeV) 


rrms (fm) 


Ps 


Pd (%) 




0.96 
0.99 


0.40 
3.22 


4.57 
2.00 


















1.00 


5.13 


1.65 


















1.10 


83.53 


0.58 


















0.80 
0.90 


3.82 
19.40 


1.86 
1.04 


96.33 
94.34 


3.67 
5.66 


1.15 
1.20 


0.77 
1.89 


3.42 
2.35 


94.89 
93.01 


5.11 
6.99 




1.00 


59.74 


0.74 


90.03 


9.97 


1.40 


12.69 


1.10 


88.10 


11.90 




1.05 


90.87 


0.66 


86.20 


13.80 


1.50 


22.91 


0.88 


86.44 


13.56 




0.80 


14.13 


1.01 


















0.90 


13.58 


1.07 


















1.00 
1.10 


34.00 
83.78 


0.77 
0.56 


















0.90 


0.56 


3.99 


99.76 


0.24 












1.00 


7.53 


1.41 


99.59 


0.41 














1.10 
1.20 


22.97 
43.80 


0.94 
0.76 


99.58 
99.58 


0.42 
0.42 













TABLE VIII: Comparison of the numerical results of the system H<,Hc in the OBE model and OPE model. 



The S^Sp(S^) systems are similar to ScSc(Sc) and the results are listed in Figs. 4 (f)-(h) and Tables VII- VIII. 
Among the six bound states, is the most interesting one. As shown in Fig. 4 (f), the r] exchange does not 

contribute to the total potential. The tt exchange is repulsive. So the dominant contributions are from the a, co, p 
and exchanges, which lead to a deep well around r = 0.6 £m and a loosely bound state. When we increase the 
cutoff from 1.50 GeV to 1.90 GeV, the binding energy of varies from 0.18 MeV to 5.58 MeV, and the RMS 

radius varies from 5.52 fm to 1.74 fm. This implies the existence of this loosely bound state. If we consider the 
TT exchange alone, only the state is bound. The percentage of the "^^i component is more than 86% when 

1.15 GeV < A < 1.50 GeV as shown in Table VIII. 



The ric^ci^c) case is quite simple. Only the rj, a and </> exchanges contribute to the total potentials. The shape of 
the potential of ^'n^o^ is similar to that of ^^'h' ■ '^^^ binding energy of this state is very small. For the spin-triplet 
SlcClc system, its S wave percentage is more than 99%. In other words, the S— D mixing effect is tiny for this system. 

We give a brief comparison of our results with those of Refs. [36, 37] in Table X. In Ref. [36], Froemel et al. 
deduced the potentials of nucleon-hyperon and hyperon-hyperon by scaling the potentials of nucleon-nucleon. With 
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A (GcV) E (McV) 



(fm) 



A (GeV) E (MeV) rrms (fm) Ps : Pd (%) 



[0,1] 



0.96 


1.07 


3.04 


0.98 


2.67 


2.08 


1.00 


4.51 


1.69 


1.20 


5.92 


1.59 


1.70 


19.88 


1.15 


0.90 


13.12 


1.06 


0.97 


4.34 


1.70 


1.00 


5.01 


1.62 


1.10 


20.96 


0.94 


1.20 


108.50 


0.48 



[0,3] 



[0,1] 



[0,3] 
S2cS2c 



0.80 
0.88 
1.00 
1.10 
1.20 



6.92 
3.05 
9.77 
26.22 
47.23 



1.53 
1.98 
1.23 
0.86 
0.72 



99.64 0.06 

99.96 0.04 

99.90 0.10 

99.79 0.21 

99.53 0.47 



TABLE IX: Numerical results of the Qc^~lc and Oc^^c systems. 



the nucleon-nucleon potentials from different models, they discussed possible molecular states such as Scc-/V, ScScc, 
EcEc etc.. The second column of Table X shows the binding energies corresponding different models while the last 
column is the relevant results of this work. One can see the results of Ref. [36] depend on models while our results 
are sensitive to the cutoff A. 



models 


Nijm93 Nijml Nijmll AV18 AV8' 


AV6' 


AV4' 


AVX' 


AV2' 


AVI' 


This work 


[ScSc]/=0 




* 


71.0 457.0 - 


0.7 


24.5 


9.5 


12.8 




1.22 ~ 21.50 


[EcEc]/=2 


66.6 




41.1 










0.7 




PcSc]/=l 




* 


53.7 




7.3 


2.8 


8.3 


0.7 


0.11 ~ 31.35 


[EcEc]/=o 


* 


* 


285.8 * 16.1 


10.8 


87.4 


53.3 


58.5 


0.7 


1.80 ~ 75.75 



TABLE X: The comparison of the binding energies of H'H' and EcEc systems in this work and those in Ref. [36]. The unit 
is MeV. "-" means there is no bound state and "*" represents exiting unrealistic deeply bindings (1 ~ 10 GeV). 



V. CONCLUSIONS 



The one boson exchange model is very successful in the description of the deuteron, which may be regarded as a 
loosely bound molecular system of the neutron and proton. It's very interesting to extend the same framework to 

investigate the possible molecular states composed of a pair of heavy baryons. With heavier mass and reduced kinetic 
energy, such a system is non-relativistic. We expect the OBE framework also works in the study of the heavy dibaryon 
system. 

On the other hand, one should be cautious when extending the OBE framework to study the heavy baryonium 
system. The difficulty lies in the lack of reliable knowledge of the short-range interaction due to the heavy baryon 
and anti-baryon annihilation. However, there may exist a loosely bound heavy baryonium state when one turns off 
the short-range interaction and considers only the long-range one-pion-exchange potential. Such a case is particularly 
interesting. This long-range OPE attraction may lead to a bump, cusp or some enhancement structure in the heavy 
baryon and anti-baryon invariant mass spectrum when they are produced in the e+e— annihilation or B decay process 
etc. 

In this work, we have discussed the possible existence of the AcAc(Ac), ScSc(Sc), ScSc(Sc), SgS^(S^) and flc^d^c) 
molecular states. We consider both the long range contributions from the pseudo-scalar meson exchanges and the 
short and medium range contributions from the vector and scalar meson exchanges. 

Within our formalism, the heavy analogue of the H dibaryon ^I'lvli '^'^^^ ^'^^ exist though its potential is attractive. 
However, the ^"1^'^ and bound states might exist. For the SgSc system, there exists a loosely bound state 

vlf[}'^l with a very small binding energy and a very large RMS radius around 5 fm. The spin-triplet state; 'J'^''3 
may also exist. Its binding energy and RMS radius vary rapidly with increasing cutoff A. The qualitative properties 
of and are similar to those of *|^'^ • They could exist but the binding energies and RMS radii are 

unfortunately very sensitive to the values of the cutoff parameter. 

For the ScSc, ScSc, S^S'^, S^S'^, ilc^^c and ilcdc systems, the tensor forces lead to the S — D wave mixing. There 
probably exist the EgEc molecules ^E'^cSc *EcEc ^^^y- ^c^c system, the cj and p exchanges are crucial 
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to form the bound states ^'^ and . If one considers the tt exchange only for the S^S^ system, there 

may exist one bound state . 

The states and are very interesting. They are similar to the deuteron. Especially, and "^-1-, 

have the same quantum numbers as deuteron. For "^-'^ , the S — D mixing is neghgible whereas for deuteron such 
an effect can make the percentage of the D wave up to 4.25% ~ 6.5% [22, 24, 38]. The D wave percentage of ^'E'S' 
is 2.12% ~ 3.17%. 

The other two states ^i-'- and ^'B/ij, arc very loosely bound S wave states. Remember that the binding energy 
of deuteron is about 2.22 MeV [39] with a RMS radius r^ms ~ 1-96 fm [40]. The binding energy and RMS radius of 
^s'S' quite close to those of the deuteron. In contrast, the state "fslH,: much more loosely bound. Its binding 
energy is only a tenth of that of deuteron. 

However, the binding mechanisms for the deuteron and the above four bound states are very different. For the 
deuteron, the attraction is from the n and vector exchanges. But for these four states, the tt exchange contribution 
is very small. Either the a (for ScSc) or vector meson (for S^S^) exchange provides enough attractions to bind the 
two heavy baryons. 

Although very difficult, it may be possible to produce the charmed dibaryons at RHIC and LHC. Once produced, 
the states ScSc and S(,S^ are stable since Sc and S', decays cither via weak or electromagnetic interaction with a 
lifetime around lO^^'^s [31]. On the other hand, Sc mainly decays into A+tt. However its width is only 2.2 MeV [31]. 
The relatively long lifetime of allows the formation of the molecular states "S"^'^ and ^'^'^ • These states may 
decay into EcA+tt or A+A+tttt if the binding energies are less than 131 MeV or 62 MeV respectively. Another very 
interesting decay mode is EccN with the decay momentum around one hundred MeV. In addition, a baryonium can 
decay into one charmonium and some light mesons. In most cases, such a decay mode may be kinetically allowed. 
These decay patterns are characteristic and useful to the future experimental search of these baryonium states. 

Up to now, many charmonium-like "XYZ" states have been observed experimentally. Some of them are close to the 
two charmed meson threshold. Moreover, Belle collaboration observed a near-threshold enhancement in e+e~ AcAc 
ISR process with the mass and width of m = (4634t^(stai.)+^(.s?/s.)) McV/c^ andTtot = {92tt'i{stat.)tl°{sys.)) MeV 
respectively [41]. BaBar collaboration also studied the correlated leading Ac Ac production [42]. Our investigation 
indicates there does exist strong attraction through the a and oj exchange in the ApAc channel, which mimics the 
correlated two-pion and threc-pion exchange to some extent. 

Recently, ALICE collaboration observed the production of nuclei and antinuclei in pp collisions at LHC [43] . A 
significant number of light nuclei and antinuclei such as (anti)deutcrons, (anti)tritons, (anti)Helium3 and possibly 
(anti)hypertritons with high statistics of over 350 M events were produced. Hopefully the heavy dibaryon and heavy 
baryon and anti-baryon pair may also be produced at LHC. The heavy baryon and anti-baryon pair may also be 
studied at other facilities such as PANDA, J-Parc and Super-B factories in the future. 
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APPENDIX 
A. The functions Ho, Hi, H2 and H^ 



The functions Hq, Hi, H2 and are defined as [17] 

Ho{K, m,r) = — (e"'"'' - e"^'') _ " ™' .-Ar 



mr " ' 2mA 

iJi(A,m,r) = -J-{e-^- -e-^-)+K ^ ~ ^ e 
mr 2m'^ 



H2{A,m,r) = 1 + _ __ - M + _ 



Ar J m m'^r^ 

rr ^ (. ^ 3 Ne"™'' / 3 3 \A2e-^'- A2-m%, , -.e'^^ 

Hs{A,m,r) = 1 + — + ^^ 1+ + 1 + Ar . 35 

\ mr m'^r'' J mr \ Ar A^r^/ mr 2m'' mr 

With Fourier transformations we have 

1 rr / A ^ 

Ho{A,m,r), 



m^ + 
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QiQj 



47r 

imPr 
An 

~12^ 



Hi{A,m,r), 



H2{A,m,r), 

{f3(A, m,r) -(5, 



(36) 



B. The coupling constants of the heavy baryons and light mesons 

In the quark model we have 

_ _ 4 TOs, _ 2 mH'„ 
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'99 



3 ruq 

2 mg^ 

3 mo 



- 1 



QaAc^c — "^Qa-qqi Qcr'Ec'S.c — "^Qcrqq 



(37) 



Because nucleons do not interact directly with the (j) meson in the quark model, we can not get g^qq in this way. 
However, using the SU{3) flavor symmetry, we have g^qq = \f2gpqq. Since gpqq is related to gpNN, all coupling 
constants of heavy charmed baryons and (j) can be expressed in terms of gpNN- 



C. The dependence of the binding energy on the cutoff 



Finally, we plot the variations of the binding energies with the cutoff. 
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